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The entanglement entropy and quantum fidelity in a hard-core-boson model with nearest- and 
next-nearest-neighbor interactions are studied numerically. By using exact diagonalization and the 
density matrix renormalization group, the eflects of interactions on entanglement entropy and fidelity 
susceptibility in the model are investigated. We focus our attention on looking for three quantum 
phase transitions. It is found that the first derivative of the entanglement entropy can indicate all of 
three phase transitions, while the fidelity susceptibility cannot predict the transition from superfluid 
to bond-order in a finite-size system. 
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I. INTRODUCTION 

In condensed matter physics, quantum phase transi- 
tions, which happen at zero temperature, have received 
much attention. In a quantum phase transition the 
properties of a many-body system will change dramati- 
cally, when a controlling parameter changes across criti- 
cal point [1] . The phase diagram in the one-dimensional 
spinless fermion model has been determinded by many 
researchers. It is well known that a transition from Lut- 
tinger liquid to charge-density-wave phase occurs when 
the repulsive interaction is twice the hopping interac- 
tion. However, long-range Coulomb interactions some- 
times also play an important role in frustrated systems 
due to their geometrical structure. For example, a spin- 
less fermion, including the nearest-neighbor and next- 
nearest-neighbor repulsive interactions can describe the 
low-energy electronic state of CuO double chains [2], 
and the repulsions plays a significant role in the charge- 
ordered-insulator-metal transition. 

The phases of spinless fermion including the nearest- 
neighbor as well as next-nearest-neighbor repulsive in- 
teractions are more complicated than one might expect. 
The model can be described by the following Hamilto- 
nian: 

H = Y.^-Mc^+l + h.c.) + Vi{cU - - \) 

i 

^V2{c\c,-]^){c\^^C,+2~\)i (1) 

where cj(ci) is the fermionic creation (annihilation) op- 
erator at site J, i is the hopping amplitude, and Vi and 
V2 are the nearest-neighbor and next-nearest-neighbor 
Coulomb interactions. The fermionic Hamiltonian ([T]) 
can be mapped onto the identical Hamiltonian for hard- 
core bosons, which is given by 

H - Y.^-t{h%+^ + h.c.) + V,{blh - l){bU,b,+i - I) 

i 

+ V2{blb,-'l)ibl2b^+2-l)], (2) 

where bl{bi) is the hardcore boson creation (annihila- 
tion) operator at site i, and t is the hopping integral and 



taken as the unit of energy in the paper. Vi and V2 are 
the nearest-neighbor and next-nearest-neighbor inter- 
actions. This mapping can be done by using a Jordan- 
Wigner transformation, as both models have the same 
spectrum. The phase diagram for hard-core bosons is 
relevant to the corresponding spinless-fermion model, as 
is shown in Refs. It is noted that the Luttinger- 

liquid phase in the spinless-fermion model is relevant to 
the superfluid phase (SF) in the hard-core-boson model, 
and other phases have the same names in the hard-core- 
boson model and the spinless-fermion model. We also 
focus our analysis on the model ^ at half filling. The 
Hamiltonian of the case V2 — 0, can be solved exactly 
via Bethe ansatz. We find a transition at Vi/t — 2 from 
the SF to the (• • • 01010101 • • • ) charge-density-wave 
(CDW-I) phase as Vi is increased. When V2 7^ 0, the 
phase structure becomes rich. It is reported that there 
exist bond-order (BO) phase[3'-'3], which is somehow 
was missed in some research ^i7j. When V2 is large, 
(• • • 00110011 • • • ) charge-density-wave phase(CDW-II) 
with a four-site unit cell appears. Here, 1(0) denotes 
the presence (absence) of a particle at a particular site. 



Recently, due to the cross over field between quan- 
tum many-body theory and quantum-information the- 
ory, the ground-state entanglement entropy and fidelity 
have been used to qualify quantum phase transitions 
in the one-dimensional spin systems [SHIB]. Our goal 
is to combine these newly developed observables with 
accurate numerical calculations, and independently de- 
termine the phase diagram of the above model. 



In the paper, we calculate the ground-state entangle- 
ment entropy and fidelity in the half-filled hardcore bo- 
son model including the nearest-neighbor as well as the 
next-nearcst-ncighbor repulsive interactions, and use 
them as good tools for searching for the phase tran- 
sition points. The remainederof the paper is organized 
as following. The measurements and the details of the 
methods to obtain the ground state are introduced in 
Sec. II. The results for entanglement entropy are pre- 
sented in Sec. HI and for fidelity in Sec IV. Finally a 
discussion is given at last in Sec. V. 
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II. MEASUREMENTS AND METHODS 

The entanglement entropy can be chosen as a mea- 
surement of the bipartite entanglement, which can be 
used to detect a quantum phase transition. The entropy 
is defined as follows. Let \g-s.) be the ground state of 
a chain system, which can be divided into two parts A 
and B. The reduced density matrix of part A can be ob- 
tained by taking the partial trace over system B, which 
is given by 

PA = TrB{\g.s.){g.s.\). (3) 
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The bipartite entanglement between parts A and B can 
be measured via the entanglement entropy as 



Sab = ~Tr{pA \og2 Pa)- 



(4) 



FIG. 1. (Color online) The entanglement entropy of two 
central qubits (a) and its first derivative (b) are plotted as 
functions of V2 for Vi = 4— V2 with t = 1 for different system 
sizes. Here and in the following graphics all quantities are 
dimensionless. 



Another concept from quantum information theory, 
ground-state fidelity, can be applied to capture the ex- 
istence of the quantum phase transitions. The general 
Hamiltonian of a quantum many-body system can be 
written as H{\) = Hq + XHj, where Hi is the driving 
Hamiltonian and A denotes its strength. If p{X) rep- 
resents a state of the system, the ground-state fidelity 
between p{\) and p{\ + 5) can be defined as 



F(A, 5) = rr[y^pi/2(A)p(A-fJ)pi/2(A)] . 



(5) 



For a pure state p = Eq- [sjcan be rewritten as 

F(A,(5) = + 5))\. F{X,S) reaches its max- 

imum value Fmax = 1 at (5 = 0. By expanding the 
fidelity in powers of S, and since the first derivative 
^^^of ~ ^' fidelity can be written as 



F{X,S) ~ 1 + 



d^F{X,S) 
2dS^ ' 



(6) 



Therefore, the average fidelity susceptibility x(A, S) can 
be given by [H] 



X(A, 5) = lim 



2[1-F(A,<5)] 



N6^ 



(7) 



where N is the length of the system. 

As is well known, it is not easy to calculate the 
ground-state fidelity and entanglement entropy, because 
of the lack of knowledge of the ground-state wave func- 
tion. The Hamiltonian ^ is not exactly solvable except 
for V2 = 0; we resort to exact diagonalization to obtain 
the ground state for small size, i.e., N — 20. For large 
size, the method of the density-matrix renormalization- 
group(DMRG) [171 [HI can be applied to obtain the 
ground state of model (2]). The total number of density 
matrix eigenstates held in a system block is m = 250 
in the basis truncation procedure. The Matlab code 
for the finite-size density-matrix renormalization group 
with double precision is used to evaluate the model with 
system sizes up to N = 100, and the truncation error 
is smaller than 10~^. Both the above two methods are 
used to calculate the model ^ with open boundary 
conditions. 



III. ENTANGLEMENT ENTROPY 

In order to avoid boundary effects, the entanglement 
entropy of two central two neighboring sites of a large 
system is calculated. By use of exact diagonalization 
and the DMRG, the entanglement entropy of the cen- 
tral two particles is plotted as a function of interaction 
for sizes N = 20, 40, 60, 80, 100 in Fig. [TJa), and the cor- 
responding first derivatives are shown in Fig. [l|b) . It is 
found that the entanglement entropy initially increases 
rapidly with increases of ¥2- When the interaction V2 
increases further, the entanglement entropy increases 
slowly, and reaches a peak value at V2 = 2.80, which 
is independent of system size. As we know from Refs. 
[21 H] , the peak of local entanglement does not coincide 
with the critical point. Considering the energy of the 
system in terms of the reduced density matrix of two 
spins at positions i and j, the energy reads 



(8) 



where Hij is the reduced Hamiltonian of the two sites 
at positions i and j, and the sum is over the total Hamil- 
tonian of the system, pij is the reduced density matrix 
of two sites at positions i and j. Moreover, the second- 
order derivative of the energy is related to the derivative 
of the reduced density matrix. This is shown by 



9^ 



^ dV2 dV2 '^'^ ' 



It is noted that a discontinuity in the second derivatives 
of the energy requires the divergence of at least one of 
the derivatives at the critical points [12] .The phase 
transition points sometimes are determined by the po- 
sitions at which the first-order derivative of the two-site 
entanglement entropy takes on its maximal and mini- 
mal values [inHHl- Here, the CDW-II to BO transi- 
tion is a second-order transition, so it should be cap- 
tured by the first-order derivative of the two-site en- 
tanglement entropy. The phase transitions of SF to 
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FIG. 2. (Color online) The scaling behavior of the extreme points of V2 versus N''. The red lines are fixed lines. 
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result obtained from the energy gap and larger than 
the result V2 —2.55 obtained from the static structure 
factor, which is also shown in Ref. [3]. 



FIG. 3. (Color online) The phase diagram of the model 
(|2| at half filling in the tilted V1-V2 plane. The dashed line 
corresponds to V2 = Vi/2, which is the boundary of the 
CDW-I and CDW-II phases when t = 0. 



BO and BO to CDW-II are infinite-order Berezinskii- 
Kosterlitz-Thouless (BKT) transitions, which can also 
be predicted by the first-order derivative of the two-site 
entanglement entropy. The first derivative of the entan- 
glement entropy is shown in Fig. [ijb). There are two 
peaks, whose locations move to high V2 with increase 
in system size, and the location of the valley moves to 
lower V2 with system size increases. These behaviors 
capture the CDW-I to SF, SF to BO, and BO to CDW- 
II phase transitions, respectively. 

The positions of the two peaks and one valley in 
the first derivative of the two-site entropy can be ex- 
trapolated to the thermodynamic limit. We investigate 
the size dependence of the positions of their extreme 
points. The changes of the critical points can be fitted 
by the formula V2 ~ V2 + aN'', where a, b are con- 
stants and N is the size of the system. We plot the 
locations of the two peaks and one valley as functions 
of N''. The fitting results are shown in Fig. [2] Here 
we obtain that = 0.77, ai = 3.012, h = -0.45, 
V^^^ = 2.60, a2 = -6, 564, 62 = -1.02, and ¥2^^ = 2.978, 
aa = -0.9133, 63 = -0.621. 

In Fig. [3j we also show the phase diagram of the 
model, which is detected via the entanglement entropy. 
It is seen that our results for the CDW-I to SF boundary 
and the CDW-II to BO boundary are consistent with 
the results in Refs. [Hill]- The difference is in the BO 
to SF boundary. The region of the BO phase is smaller 
than the region captured by the static structure factor. 
For Vi — 4, we capture the BO to SF critical phase 
transition point V2 — 2.90. It is consistent with the 



IV. FIDELITY 

We apply the DMRG to calculate the ground-state 
fidelity susceptibility with system size N up to 100 and 
5 = 0.01. The ground-state fidelity susceptibility x is 
plotted as a function of the interaction V2 for different 
sizes in Fig. |4]in dimensionless units. Two peaks in the 
fidelity susceptibility is found. The peaks of x increase 
when the system size increases. One of the peaks' loca- 
tions moves to high V2 up to a particular value as the 
system size increases, and the other moves to a slightly 
lower V2 down to a particular value as the system size 
increases. The fidelity measures the similarity between 
two states, while quantum phase transitions are intu- 
itively accompanied by an abrupt change in the struc- 
ture of the ground-state wave function. This primary 
observation motivates researchers to use the fidelity to 
predict quantum phase transitions. 

A peak in the fidelity susceptibility indicates a quan- 
tum phase transition. Here, scaling of the extreme 
points of the fidelity susceptibility as the system length 
increases is also investigated. We find that the changes 
of the maximal points can also be fitted by the for- 
mula V2 ^ V2 + aN^, where a,b are constants and N 
is the number of the system. The results for the loca- 
tions of the fidelity susceptibility can be used to inves- 
tigate the quantum phase transition points in the ther- 
modynamic limit. We plot the locations of the max- 
imum fidelity susceptibility as a function of and 
show the numerical fit in Fig. [4] (b) and |4] (c) . We 
obtain that = 0.76, ai = -1.522, 61 = -0.262, and 
¥2^^ = 2.978, 03 = -0.9133, 63 = -0.621. 

The quantum phase transitions points are partially 
consistent with the results obtained from the entangle- 
ment entropy. It is shown that the fidelity susceptibility 
can predict the CDW-I to SF transition and the BO to 
CDW-II transition, however, it cannot predict SF to BO 
transition. As we know, the BO to CDW-II and SF to 
BO transitions are both infinite-order BKT transitions. 
The Numerical results clearly show that the BKT tran- 
sition point is detectable to the BO to SF transition 
by divergence of the fidelity susceptibility [22,. Our 
results are similar to the results, which the entangle- 
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FIG. 4. (Color online) (a) Fidelity susceptibility versus 
interaction V2 for Vi = 4 — V2 with t = 1 for different system 
sizes. (b),(c) The scaling behavior of the extreme points of 
V2 versus A''*'. The red lines are fixed lines. 



perturbations, where the fidelity susceptibility may not 
diverge. It would be interesting to use infinite-time- 
evolving block decimation[27^ to calculate the ground- 
state fidelity, and use it to detect the superfluid to bond 
order transition. Here we do not show the phase dia- 
gram of the model which can be detected by using the 
fidelity susceptibility. 



V. DISCUSSION 

In the paper, we employ the entanglement entropy 
to study numerically the hard-core-boson model with 
nearest and next-nearest-neighbor interactions. By us- 
ing exact diagonalization and the density-matrix renor- 
malization group, the effects of interactions on the en- 
tanglement entropy in the model are investigated. It 
is found that the first derivative of the entanglement 
entropy can detect all the quantum phase transitions 
in the model, and the phase diagram is given. We also 
test the ability of the fidelity susceptibility to detect the 
phase transitions. It is found that the fidelity suscep- 
tibility can predict charge-density-wave II to bond or- 
der and charge-density-wave I to superfluid transitions, 
however, the fidelity susceptibility cannot predict the 
superfluid to bond order transition for finite-size sys- 
tems. 



ment entropy predicts the BO phase in the extended 
Hubbard model [24], and the fidelity can not [25]. The 
BKT phase transition point sometimes cannot be well 
characterized by the ground-state fidelity for finite-size 
systems [53]. A other reason may be that, according 
to the analysis of Ref. [26] , a divergence in the fidelity 
susceptibility implies a quantum phase transition, but 
the converse is not true. This means that there are 
quantum phase transitions driven by particularly weak 
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